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Magnetic  Shear  Stabilization  of 
Diocotron  Instability 

S.  Kondoh,  T.  Tatsuno,  and  Z.  Yoshida 


Graduate  School  of  Frontier  Sciences,  University  of  Tokyo, 
7-3-1  Hongo,  Bunkyo-ku,  Tokyo  113-0033,  Japan 


Abstract.  The  diocotron  instability  in  a magnetized  non-neutral  plasma  is 
a close  cousin  of  the  Kelvin-Helmholtz  instability.  A sheared  magnetic  field 
brings  about  coupling  between  the  diocotron  modes  and  the  Langmuir  waves  that 
propagate  along  the  magnetic  field.  Motion  of  electrons  parallel  to  the  magnetic 
field  cancels  the  electric  charge  produced  by  the  diocotron  modes,  resulting  in 
stabilization  of  the  diocotron  instability. 

I INTRODUCTION 

Recently  a variety  of  new  concepts  on  non-neutral  plasma  confinement 
has  been  proposed  [1-3],  which  significantly  differ  from  the  conventional  Pen- 
ning/Malmberg  trap  [4].  The  Prototype  Ring  Trap  (Proto-RT)  experiment 
[1,5]  is  aimed  at  pure  magnetic  confinement  of  a toroidal  non-neutral  plasma 
that  is  not  in  a rigid-rotating  thermal  equilibrium  state.  In  such  a system, 
the  plasma  flow  is  generally  sheared,  and  the  diocotron  instability  [6]  can  be 
destabilized.  Application  of  magnetic  shear  is  expected  to  be  most  effective 
to  stabilize  the  electrostatic  modes.  However,  exact  stability  analysis  has  not 
been  completed,  except  for  the  special  case  of  an  electron  beam  with  a rela- 
tivistic speed  [7]. 

The  physical  mechanism  of  the  diocotron  instability  is  explained  as  follows 
[8]  (see  Fig.  1):  When  a non-neutral  slab  plasma  has  a finite  thickness,  a 
perturbation  on  one  of  the  two  plasma  surfaces  produces  surface  charges.  The 
resulting  perturbed  electric  field  yields  an  E x B flow  in  the  plasma,  and  the 
opposite  surface  is  also  perturbed.  The  motion  of  the  opposite  surface  brings 
about  a reciprocal  perturbation,  and  the  waves  on  the  two  surfaces  couple  with 
each  other.  Under  certain  conditions,  this  coupling  yields  a positive  feedback, 
and  the  diocotron  instability  occurs. 

When  a non-neutral  plasma  is  confined  in  a uniform  magnetic  field,  the 
diocotron  modes  propagating  in  the  perpendicular  direction  to  the  magnetic 
field  are  independent  of  any  modes  that  propagate  in  the  parallel  direction. 
However,  if  the  magnetic  field  has  a shear  (see  Fig.  2),  the  wave  vector  may 
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FIGURE  1.  Physical  picture  of  diocotron  modes  in  a uniform  magnetic  field.  Perturbation 
on  one  of  the  two  plasma  surfaces  produces  the  surface  charge  and  causes  the  electrostatic 
field  perturbation.  This  perturbed  electric  field  shakes  the  body  of  the  plasma  itself  through 
Ex  B drift,  and  the  other  surface  is  also  perturbed.  The  perturbation  on  the  latter  surface 
in  turn  shakes  the  former  one  in  the  same  way.  Thus,  the  waves  on  the  two  surfaces  couple 
with  each  other.  Under  certain  conditions,  the  diocotron  modes  can  be  unstable. 
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FIGURE  2.  Physical  picture  of  stabilizing  effect  of  a sheared  magnetic  field  on  diocotron 
modes.  The  wave  vector  almost  always  has  a local  parallel  component,  and  the  diocotron 
modes  cannot  be  independent  of  the  parallel  modes,  such  as  the  Langmuir  wave  or  the 
plasma  oscillation.  Therefore,  a coupling  between  them  is  caused.  In  a cold  non-neutral 
plasma,  the  surface  charge  perturbation  produced  by  the  diocotron  modes  is  canceled  by 
the  parallel  motion  of  charged  particles.  Thus,  the  diocotron  instability  is  stabilized  by  the 
sheared  magnetic  field. 


have  a local  parallel  component  /i:||(a:),  and  the  diocotron  modes  interact  with 
the  parallel  modes,  such  as  the  Langmuir  wave  or  the  plasma  oscillation.  In 
a cold  non-neutral  plasma,  the  surface  charge  perturbation  produced  by  the 
diocotron  modes  is  short-circuited  by  the  parallel  motion  of  charged  particles, 
if  the  diocotron  frequency  u)d  is  much  smaller  than  the  plasma  frequency 
Wp,  i.e.,  when  a low  density  plasma  embedded  in  a strong  magnetic  field. 
Therefore,  we  expect  that  the  diocotron  instability  is  stabilized  in  a sheared 
magnetic  field. 

In  this  paper,  we  consider  a slab  plasma  with  a flat-top  density  profile 
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and  show  the  stabilizing  effect  of  nnagnetic  shear  analytically.  The  diocotron 
instability  is  formally  equivalent  to  the  Kelvin-Helmholtz  instabilities  in  fluids 
and  plasmas  [9].  The  magnetic  shear  stabilization  of  these  instabilities  is  of 
the  common  interest  and  has  a variety  of  applications  (see  Sec.  IV). 

II  EIGENEQUATION  FOR  DIOCOTRON  MODES 
IN  A SHEARED  MAGNETIC  FIELD 

A Slab  plasma  model  in  a sheared  magnetic  field 

We  consider  a slab  electron  plasma  embedded  in  a sheared  magnetic  field 
(see  Fig.  2).  The  plasma  has  a finite  thickness  2A  in  the  x-direction.  We 
assume  that  all  equilibrium  quantities  are  functions  of  only  x.  We  consider  a 
sheared  magnetic  field  such  as 

B = (1) 

where  B^  is  a constant.  Since  a non-neutral  plasma  has  a self-electric  field, 
there  is  a stationary  flow  that  is  approximately  equivalent  to  the  E x B drift 
for  low  densities. 

The  governing  equations  are 

dti 

— + u • Vn  + nV  • V = 0,  (2) 

^ + („.V)u  = -l(^-br  xB),  (3) 

= (4) 

where  n is  normalized  by  the  typical  electron  density  no(0),  t by  the  inverse  of 

the  diocotron  frequency  — £oBz/eno(0)  (eo  is  the  vacuum  dielectric  con- 
stant and  e is  the  elementary  electric  charge),  the  spatial  coordinates  a;,  y,  x 
by  the  half  thickness  of  the  slab  plasma  A,  v by  the  mean  flow  velocity  at 
the  plasma  surface  |wo(l)|,  B by  the  axial  magnetic  field  E by  the  mean 
electric  field  at  the  plasma  surface  |£'(1)|,  s = Wp/wc  = 3,  dimen- 

sionless parameter,  and  E = —'V(f>.  Since  we  consider  a low-density  plasma 
in  a strong  axial  magnetic  field,  we  may  assume  s <C  1.  In  this  limit,  we  can 
replace  Eq.  (3)  with 


-Vd>  X B{x) 
B{x)^  ’ 


9u||  , , 1 

-^  + (v  • V)U||  = -V||d-, 


(5) 

(6) 
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where  B{x)  = By(xY  + Linearizing  Eqs.  (2),  (5),  and  (6),  we  finally 
obtain  [10] 

( ,0,  \ 1 d fkxno\  . 


fc  (/>i  I + 


uj  — fcitJxo  dx  \ B 


(f>i  + r C^) 

s^(u;  - AcxUio)^ 


where  we  have  Fourier-transformed  all  perturbed  variables  y,  z,  <)  as 

'H{x,y,z,t)  = ^'o(a:}  -h  'l'i(a-)  exp[/(u;/  - k^y  - k^z)].  {i 

III  STABILIZING  EFFECT  DUE  TO  PARALLEL 

MOTION 


A Non-resonant  frequency  regime 

First,  we  show  that  the  diocotron  modes  are  stabilized  (o;j  = 0)  for  wave 
numbers  without  resonance  between  the  phase  velocity  and  the  plasma  flow, 
that  is,  Wf  — kj_vj_o  ^ 0 for  all  x.  Multiplying  Fq.  (7)  by  <j>*  and  integrating  it 
over  (—00,00),  we  obtain  from  the  imaginary  part 


/■°°  r 1 d 
' J-00  \u>  — /cxt’xoj  dx 


1 d / A;xno\  2noA:j|(a;r  — ^x'^^xo) 

k±Vio\dx  \ B ) s^ltu  — A:xUxo|‘* 


|(/)i \'^dx  = 0. 


Here  we  used  the  boundary  condition 


<;^i(±oo)  = 0.  (10) 

Since  s^  1 and  tUp  — fcx^’xo  0 at  any  point  in  the  plasma  region,  we  obtain 
wj  = 0,  which  means  stability.  This  mathematical  treatment  is  same  as  the 
standard  Rayleigh’s  analysis  [11]. 


B Dispersion  relation  with  resonances 

If  the  plasma  has  a resonant  point,  the  analysis  in  Sec.  Ill  A does  not  apply 
to  check  whether  the  eigenvalues  cu  for  Eq.  (7)  are  real  or  not.  In  this  case, 
we  need  to  solve  Eq.  (7)  directly.  The  eigenfunction  determined  by  Eq.  (7) 
is  oscillatory,  because  the  sign  of  the  last  term,  which  we  assumed  to  be  very 
large  (s  <g  1),  is  positive.  If  uj  is  not  real,  the  real  and  imaginary  parts  of  the 
eigenfunction  have  a relative  phase  angle  of  about  7t/2.  When  we  consider  a 
density  profile  with  a sharp  boundary,  we  have  to  connect  both  the  real  and 
imaginary  parts  of  the  eigenfunction  at  the  plasma  surfaces  using  the  same 
boundary  condiotion.  If  both  of  them  have  a different  phase  angle,  this  process 
fails,  which  implies  that  u;  must  be  real. 


.^14 


The  essential  characteristic  of  this  eigenvalue  problem  is  well  understood  by 
the  following  simplified  model.  First,  we  neglect  the  second  term  in  the 
bracket  of  the  first  term  of  Eq.  (7)  in  the  plasma  region  ( — 1, 1),  since  it  is 
much  smaller  than  the  last  term  when  no  ~ 1.  We  also  assume  that  k^UolB 
jumps  at  a;  = ±1  and  its  variation  is  negligible  anywhere  else,  i.e. 

+ (11) 


where  f{x)  is  a given  finite  function.  Furthermore,  we  assume 
— = const.  > 1, 


(12) 


kx(a:)uxo(2:)  = x. 


(13) 


Under  these  assumptions,  we  can  solve  Eq.  (7)  analytically,  and  the  dispersion 
relation  is  given  by  [10] 

(«  + kuj-.y  - eul  +(^k+~y  + /(-!)/(!)  -{k  + D (/(-I)  + /(I)) 


2ai 


+ku)v{f{-l)  - f{l))  + i{a  + fco;i)(/(-l)  - /(I)  - 2fcu;r)  (w  + 1) 


(a  - kco-,y  - kW,  + U + - + /(-1)/(1)  - U + r (/(-I)  + /(!)) 


+ka;r(/(-l)  - /(I))  - i{a  - ku>y{f{-l)  - /(I)  - 2ka)r)J(w  - 1)^“’,  (14) 

where  u>  = Wr  + iwi-  We  can  show  that  w is  real  for  Eq.  (15).  where  Wf  = Re  w 
and  o)j  = Imw.  Taking  the  absolute  number  of  Eq.  (15)  gives 


Mi 

1^2 


exp 


a arg 


c-^)] 


(15) 


where 


= (a  + ku-,y  - kW,  + {k  + 1)'  + /(-!)/(!)  -{k  + D (/(-I)  + /(!)) 

+W(/(-l)  - /(I))  + i{a  + kc,-, ){/{-!)  - /(I)  - 2ku;r),  (16) 

A,  = {a^  ku:-,y  - kW,  + (fc  + 1)'  + /(-!)/(!)  - (^  + ^)  (/(-I)  + /(!)) 

+W(/(-l)  - /(!))  - i{a  - ku-M-l)  -/(!)-  2fcu;r).  (17) 

If  Wj  > 0,  the  left-hand  side  of  Eq.  (16)  is  greater  than  unity,  while  the  right- 
hand  side  is  less  than  unity.  Therefore,  wj  > 0 cannot  be  satisfied.  If  o;j  < 0, 
the  left-hand  side  of  Eq.  (16)  is  less  than  unity,  while  the  right-hand  side  is 
greater  than  unity.  Therefore,  wj  < 0 cannot  be  satisfied.  Thus  cjj  = 0,  which 
means  stability.  If  wj  ^ 0,  the  eigenfunctions  of  the  three  regions  (jij,  4>n,  and 
cannot  be  connected  properly  at  x = ±1. 
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IV  SUMMARY 


We  have  shown  that  the  magnetic  shear  has  a strong  stabilizing  effect  on 
the  diocotron  instability.  The  fluid  motion  parallel  to  the  magnetic  field  short- 
circuits  the  charge  perturbation  of  the  diocotron  modes.  The  scaling  param- 
eter is  s = cVQ/tJp.  Since  the  time  .scale  of  the  parallel  motion  is  ~ Wp',  the 
condition  5 1 enables  the  parallel  motion  of  the  plasma  to  cancel  the  per- 

turbed charge  sufficiently.  Typical  non-neutral  plasmas  in  laboratories  satisfy 
this  condition. 

Mathematically  the  last  term  of  the  eigenequation  (7)  prohibits  non-real 
eigenvalues,  because  the  last  term  makes  the  eigenfunction  oscillatory.  If 
oj  ^ TZ,  the  relative  phase  angle  between  the  real  and  imaginary  parts  of 
the  eigenfunction  is  about  7t/2.  This  phase  angle  disables  both  the  real  and 
imaginary  parts  of  the  eigenfunction  to  be  connected  properly  at  the  plasma 
surfaces. 

We  note  that  our  analysis  is  based  on  a modal  approach  which,  however, 
may  not  be  complete  for  non-Hermitian  systems  [12-14].  There  remains  a 
possibility  of  secular  algebraic  behavior,  although  we  have  shown  that  there  are 
no  exponentially  unstable  modes.  This  problem  will  be  discussed  elsewhere. 
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